Weak Charge Quantization as an Instanton of Interacting cr— model 
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Coulomb blockade in a quantum dot attached to a diffusive conductor is considered in the framework 
of the non-linear cr-model. It is shown that the weak charge quantization on the dot is associated 
with instanton configurations of the Q-field in the conductor. The instantons have a finite action 
and are replica non-symmetric. It is argued that such instantons may play a role in the transition 
regime to the interacting insulator. 
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The physics of interacting electronic systems in the 
presence of disorder has been a subject of an intense 
study already for a few decades 0]. Various theoretical 
approaches have been developed for description of both 
metallic and insulator phases. The non linear cr-model 
(NLcrM) in the replica or dynamic formulation 
has proven to be the most powerful tool to deal with the 
weakly disordered (metallic) phase. Despite of a lot of ef- 
forts, the NLcrM was not really successful in treating the 
insulating and transition regimes. On the other hand, 
approaches based on the idea of hopping conductivity ||] 
are quite efficient in the phenomenological description of 
the insulating phase. This is not the case in the theory of 
non-interacting localization, where the NLcrM provides 
a unified treatment j7| of both metallic and insulating 
regimes (at least in low dimensions, d < 2 + e). It is also 
known to give some insight into the physics of the 
integer quantum Hall transitions. 

To describe the insulating phase one usually assumes 
that the system may be divided into weakly connected is- 
lands whose charge is essentially quantized. The charge 
may be changed only by an addition or subtraction of 
an integer amount of e, which costs a finite Coulomb en- 
ergy. As a result, the important charge configurations are 
quantized and inhomogeneous. Approaching the metal- 
lic phase, the coupling between islands becomes better 
and charge quantization (Coulomb blockade) effects are 
less pronounced. Finally, deep in the metallic phase 
the charge is continuous and practically homogeneous 
(due to screening). In a perturbative treatment of NLcrM 
(which is usually involved in description of the metallic 
phase) there is no room for the discreetness of the elec- 
tron charge. It is doubtful, therefore, that any NLcrM 
perturbative approach can provide a picture of the tran- 
sition or insulating regimes, where charge quantization is 
of primary importance. The purpose of this Letter is to 
show how the weak residual charge quantization may be 
incorporated in the diffusive NLcrM formalism. 

To this end I consider a toy problem of a large quan- 
tum dot with the capacitive interaction connected to a 
bulk reservoir by a diffusive non-interacting wire. The 
dot is open in a sense that the wire conductance is large: 



G ^ Gq = e^/{2-KK). One may, nevertheless, ask 
whether any Coulomb blockade effects can be seen on 
such a dot. In particular, we shall look for the oscilla- 
tory part of the electron number, N{q), on the dot as a 
function of the continuous background charge, q. (For a 
closed dot, G ^ 1, such oscillations are, of course, ex- 
tremely pronounced ||^.) In a limiting case of a short (see 
below) wire this problem was recently solved by Nazarov 
[10|. Employing the model introduced earlier by Matveev 
he derived an expression for an amplitude, A, of the 
Coulomb blockade oscillations for a dot connected to a 
reservoir by an arbitrary (delay free) scattering region: 
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Ri 
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where Ri are eigenvalues of the scattering region reflec- 
tion matrix. Then, e.g. N{q) = Acos(27rg -f- ry), where 
?7 is a phase, which depends on the scattering matrix 
of the contact. For a diffusive wire the eigenvalues Ri 
are random functions of the disorder realization with the 
probability density P{R) = G/[2(l - R)^/R], where 
G — J dR{l — R)P{R) is the average conductance of the 
wire (hereafter conductance is measured in units of Gq). 
Employing Eq. (|l]), one finds for the typical oscillation 
amplitude 

g(lnA)^g-4G^ (2) 

where the angular brackets stand for disorder averag- 
ing. This result holds for a metallic wire with L ^ ^, 
where L is the wire length and ^ is localization length, 
and therefore G ^ 1. Moreover, to treat the wire 
as an instantaneous scattering region one assumes that 
L <t: Lt = ^riD/(2iiT) fill, where D is diffusion con- 
stant of the wire and T is temperature. Since, according 
to Eq. (||), the Coulomb blockade is exponentially sup- 
pressed, we refer to it as the weak charge quantization. 

We are interested in the wire rather than in the dot. 
One would like to show, thus, how the weak charge quan- 
tization, Eqs. (P, (||), may be understood (and extended 
for L > Lt) from the point of view of the NLcrM of the 
wire. Since the coupling constant of the NLcrM is 1/G, 
the expected result, Eq. is the non-analytic func- 
tion of the coupling constant. Therefore, it can not be 
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obtained in any perturbative treatment of NLcrM. In- 
deed, we shall show that these are the instanton con- 
figurations with the finite action, which are responsible 
for the weak charge quantization. In our toy model the 
electron-electron interactions and disorder are spatially 
separated on the dot and the wire correspondingly. Nev- 
ertheless, one has to consider the full interacting NLcrM 
of Finkel'stein with the Q-field having a structure 
both in the energy and replica spaces. The finite ac- 
tion instantons, mentioned above, involve rotations be- 
tween two or more different replicas. Thus each instan- 
ton configuration is replica non-symmetric. Eventually 
the replica symmetry is restored by instantons with all 
possible replica permutations. The model provides the 
first example of a non-perturbative and non replica- 
diagonal application of the interacting NLcrM. We con- 
firm Nazarov's result, Eq. (^), and extend it to the case 
of a long, L > Lt, wire and 2D diffusive conductor. 

Let us consider a large |jl^ quantum dot connected 
to a reservoir by a non-interacting diffusive wire. The 
Coulomb interaction on the dot has the form 



(3) 



where Ec = e^/(2C) is the capacitive charging energy, 
N is the electron number operator on the dot and q is 
the background charge. Our aim is to calculate the os- 
cillatory part of the average (over quantum fluctuations) 
number of electrons on the dot, N{q)- The next step is to 
average the result over realizations of disorder in the dif- 
fusive wire. The random phase of oscillations results in 
{N{q)) — 0. Therefore, to calculate a typical oscillation 
amplitude one has to look for the correlation function 
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{N{q)N{q')) = lim -9,9,, (Z"(<7)Z"(g')) , 



n— »o n 



(4) 



where we introduced the partition function of the dot- 
wire system, Z{q), along with the replica trick to facili- 
tate the disorder averaging. The partition function of the 
open dot may be written in a standard way as the 
imaginary time functional integral over the scalar field, 
$(r). If T > A [?], the $ field configurations may be 
classified by the integer winding number, VF, by imposing 
the boundary condition <I>(r + P) = $(r) -I- 2ttW, where 
/3 = 1/T. In our case both the field and the correspond- 
ing winding number are 2n-component replica-vectors, 
* = {$1, . . . $2n} and W = {Wi, . . . W2n}- The correla- 
tion function takes the form 

{Z^{q)r\q'))= Y^e^^h^liqY^Wa + q'Y. Wa]\ 



W=-oo 



a— 1 a— n+1 



X IV^ cxp < — 



^2n ] 



(5) 



Here Sw{^) is the average effective action which is a re- 
sult of tracing out degrees of freedom of the wire. The 
latter are connected to the dot's variable, through 
the boundary condition at the dot-wire interface (for 
T > A). In terms of the NLcrM the effective ac- 
tion may be written as 



-S„(*)^yPQ exp|^yda;Tr{-i?(VQ)2 + 4eQ}| 



(6) 



where v is the density of states (DOS) in the wire and 
£ = S'^^'SmSmk IS a matrix diagonal in replica and Matsub- 
ara spaces, e,„ — 27rT(m + i). The Q-field is also a ma- 
trix, Qg^ej. (x), both in the replica space, a,b — 1, . . .2n 
and the Matsubara space, m,k — 0, ±1, ±2, . . . 

On the wire-reservoir and the wire-dot interfaces the 
Q-field obeys the following boundary conditions: 



Q{x ^L)^ A = 5"''sign(e„)(5„fe ; 
Q(a; = 0) = e-**(^)A^_^,e'*(^') , 



(7a) 
(7b) 



correspondingly. In the last equation we employed imag- 
inary time representation of the Q-field instead of the 
frequency representation. In addition, the normalization 
condition Q^ = 1 must be obeyed at every point in space. 

For a good metal, G = hvD/L ^ 1, the Q-integral 
may be calculated in the saddle point approximation. 
The extremal Q configurations are given by the Usadel 
equation 



V(L'QVQ) - [e,Q] ==0, 



(8) 



where Q satisfies Q^ = 1 and the boundary condi- 
tions (^). Our strategy, thus, is to find the solution of 
Eq. (||) for a fixed $ and identify 5u,($) with the action 
on the extremal Q configuration. 

The oscillatory (in q ^ q') component of the correla- 
tion function originates from non-zero winding numbers, 
cf. Eq. (^). Before proceeding in this direction, we shall 
make a few remarks about the W — sector. In this 
case the functional integral on the r.h.s. of Eq. (^, may 
be minimized by the replica diagonal anzatz M 



Q^^,(a;) = e-'K(^'^)A^_^,e*^(^''^' , 



(9) 



where K ~ Ka (t, x) is replica and imaginary time di- 
agonal matrix with K(r, L) = and K(t, 0) = ^{t)- 
Substituting the anzatz (^) into Eq. (^) and minimizing 
the action with respect to K |jl^, one finds 



5„($)^^V^^^coth. 



47r/3 



E, 



1/2 



Et 



where Et = fiD / is the Thouless energy of the 
wire. In the short wire limit, |w,„| ^ Et^ one recov- 
ers the phenomenological action of an ohmic environ- 
ment |l|,|l§ Su, = G/(W) ^„ Although 
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configurations with W = do not affect the thermody- 
namics of the dot, cf. Eq. (||), they change its tunnehng 
DOS ~ (exp{i(*(T) — *(r')}).^. In the small tempera- 
ture (or short wire) limit, T <^ Et, one finds the well- 
known [ |19[ power— law DOS - T^/G. In the other limit- 
ing case, T Et, the DOS is given by exp{ — ^y^^}. 
This is nothing but the ID Altshuler-Aronov zero-bias 
anomaly ||2^,^,Q for the case of interacting dot and non- 
interacting wire. One can easily include long-range in- 
teractions inside the wire, which lead to the dynamical 
screening term in the K action , to obtain the stan- 
dard ID result ||2^. Therefore, both the phenomenolog- 
ical theory of an ohmic bath and the zero-bias anomaly 
are consequences of the W = 0, replica-diagonal sector 
of the NLcrM. 

We turn now to our main subject — the non-trivial 
winding numbers, W ^ 0. One may try to solve Eq. (j^) 
by the replica diagonal anzatz, Eq. (|^). Then the phase 
Ka{T,x) acquires Wq vortices in the {x,t) plane. As a 
result, the corresponding action is logarithmically large, 
~ GlnL/A, where A is a microscopic cutoff length. 
(Note that this logarithmic divergence continues well in 
the ballistic regime, therefore A is likely to be of the or- 
der of the Fermi wavelength.) On the other hand, one 
may find solutions of Eq. (|^) with much smaller action, 
Sw ~ G. For this one must allow rotations between dif- 
ferent replicas and impose the following condition on the 
possible W configurations: 



^ W^a = . 



(11) 



Equations (Q), and ( |lT| ) immediately lead to the fol- 
lowing conclusions: (i) {N{q)) = 0; (ii) the correlation 
function {N{q)N{q')) is a function oi q — q' only. 

To demonstrate how such solutions may be con- 
structed, let us consider the optimal realization of the 
*-field, *(r) ^ 2ttWt/(3. Employing Eq. (^), one finds 



(12) 



This is equivalent to the local, replica dependent shift of 
the chemical potential on Wa Matsubara units. Consider 
e.g. the simplest non-trivial winding number configura- 
tion consistent with Eq. ( pT| ) having only two non-zero 
components of W: 

Wa = l; a€[l,n]; Wb ^ -1; b e [n + l,2n] (13) 

and Wc = for c ^ a, 6. In this case one may satisfy the 
boundary conditions Eqs. ([7a|), ( p^ ) along with = 1 
by choosing 



cos9{x) sm9{x) 
sm9{x) —cos9{x] 



(14) 



where 9(0) = tt, 9{L) = and all other elements of Q are 
equal to those of A. The Usadel equation takes the form 
of the anharmonic pendulum equation 



Ll'^^9-sm9 = 0. 



(15) 



Its solution may be written in terms of the elliptic inte- 
grals. The corresponding action is equal to 



Sin — G- 



L 
2L^ 



-.E{k) 



1- k2L 
K Lt 



+ 0{n), (16) 



where y/KK{K) — L/Lt and K{k) and E{k) are the com- 
plete elliptic integrals of the first and second kind. In the 
two limiting cases one obtains (in the n — > limit) 



^G 



1 



4T 

ttEt 



2G{Lt), 



T-^Et\ 
T-^Et; 



(17) 



where G{Lt) = e^vD/Lr = G^J2ttT/Et is the dimen- 
sionless conductance of the length Lt of the wire. In the 
zero temperature case the instanton spreads uniformly 
over the wire, 9{x) = it{L — x)/L. In the opposite limit, 
L ^ Lt, it is confined to the region of the size Lt near 
the dot: 9{x) = 27r — 4arctan[exp{a;/LT}]. The low tem- 
perature (short wire) limit may be directly compared to 
Nazarov's result, Eq. (||) (we calculate {N'^) rather than 
(|iV|) thus the factor of two difference). In general case 
Eq. ( p^ provides the entire temperature and wire-length 
dependence of the oscillation amplitude (with the expo- 
nential accuracy). Winding number configurations other 
than Eq. ( p^ are exponentially less probable. The ana- 
lytical continuation, n ^ 0, is straightforward if one no- 
tices that there are precisely different configurations 
like Eq. (|l|) with the same action and q — q' dependence. 
Thus the factor in Eq. (^ is cancelled and the re- 
mained expression is n-independent. It is important to 
notice that the rotational symmetry between the replicas 
is broken down to the permutation one by the choice of 
the winding number configuration. As a result there is 
no degenerate saddle point manifold, which n-dependent 
volume plays a central role in the case of non-interacting 
level statistics 

One can define the current operator on the instanton 
trajectory as J = iQVQ = cry\/9(x), where <7y is the 
Pauli matrix in the space defined in Eq. (Q. In the limit 
of the long wire, L 3> Lt, the current is non-uniform 
along the wire since, V0 ^ const. By virtue of the con- 
tinuity equation, this seemingly leads to the charge ac- 
cumulation in the region x ~ Lt- Where this indeed 
the case, our model of the non-interacting wire would 
loose its validity and finite temperature results should be 
modified. This is not the case, since the current, J, is 
pure replica off-diagonal (because of the ay matrix) and 
therefore completely decouples from the Coulomb inter- 
actions. Interpretation of this counter-intuitive phenom- 
ena is that the replica structure of the instanton takes 
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care of the mesoscopic fluctuations effect. It describes 
sample specific energy (and thus temperature) depen- 
dence of the transmission coefficients, introducing retar- 
dation neglected in Nazarov approach. This retardation 
is not associated with the accumulation of any real charge 
in some portion of the wire. 

Finally, we address the 2D setup with the dot of the 
radius d placed at the center of a 2D disordered disk of 
the radius L. In this case the elementary instanton is 
given by Eqs. (|l|), (0) with 9 = e{r) which is a solution 
of the following equation 







(18) 



with 9{d) = TT and 0{L) = 0; r is the polar radius. For 
L ^ Lt the equation is solved by 



e{r) = TT 



1 - 



Inr/d 



InL/d 



with the corresponding action 



\TLL/d 



(19) 



(20) 



At higher temperature, d ^ Lt < L, with the loga- 
rithmic accuracy one must substitute Lt instead of L 
in Eq. (pO|). By analogy with Eq. ( p7| ) the final result 
may be written as — tt^G{Lt)/4:, where G{Lt) = 
2'Ke'^vD / \uLt / d is the conductance of 2D disk with the 
outer radius Lt and the inner one d. Note that the tem- 
perature dependence of the 2D action is [In _D/(c?^r)]^^ 
much slower than y/T in one dimension. 

In conclusion, we have shown that the residual charge 
quantization in the diffusive interacting systems may be 
described by the instantons of the interacting NLcrM. 
These instantons are replica non-symmetric finite action 
configurations of the Q-field, which extremize the action 
with the "twisted" boundary conditions. As an illustra- 
tion, we have calculated the (exponentially small) am- 
plitude of the Coulomb blockade oscillation in the open 
quantum dot connected to ID (2D) diffusive conductor 
of arbitrary length (radius) . Although the instantons are 
exponentially rare in the diffusive regime, they are re- 
sponsible for the charge quantization and inhomogeneity 
- the features which dominate the physics of the insu- 
lating phase. One may expect them to be increasingly 
important as the insulator regime is approached. 
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